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We investigate the wetting properties of surfaces patterned with fine elastic hairs, with an emphasis on iden- 
tifying superhydrophobic states on hydrophilic hairs. We formulate a two dimensional model of a large drop 
in contact with a row of equispaced elastic hairs and, by minimising the free energy of the model, identify the 
stable and metastable states. In particular we concentrate on partially suspended states, where the hairs bend 
to support the drop - singlet states where all hairs bend in the same direction, and doublet states where neigh- 
bouring hairs bend in opposite directions - and find the limits of stability of these configurations in terms of 
material contact angle, hair flexibility, and system geometry. The drop can remain suspended in a singlet state 
at hydrophilic contact angles, but doublets exist only when the hairs are hydrophobic. The system is more likely 
to evolve into a singlet state if the hairs are inclined at the root. We discuss how, under limited circumstances, 
the results can be modified to describe an array of hairs in three dimensions. We find that now both singlets 
and doublets can exhibit superhydrophobic behaviour on hydrophilic hairs. We discuss the limitations of our 
approach and the directions for future work. 

I. INTRODUCTION 

There is widening interest in understanding how fluid streams and drops interact with micropatterned surfaces. For example 
chemical patterning can be used to direct streams of fluid, and hydrophobic surfaces which are rough on the micron or 
nanometre length scale can exhibit superhydrophobicity, characterised by contact angles near 180° and easy roll-off 0]. A 
primary motivation to study micropatterned substrates comes from their potential applications as, for example, more efficient 
water and dew repellent materials, dehumidifiers or low drag surfaces. Inspiration for these developments has come in part from 
nature where plants, such as the lotus 0] and nasturtium, have superhydrophobic leaves, and water- walking insects, like the 
water strider yfl, have superhydrophobic legs to allow them to interact with the the water in their environment. 

Many of the natural surfaces that show strong water repellency are covered with tiny hairs. These can be long and flexible 
as on the leaves of Lady's Mantle or form a dense, spiky array as on the legs of the water strider. There have also been recent 
advances in the microfabrication of hairy surfaces 14], [3]. To fully exploit these possibilities it is important to gain a better 
theoretical understanding of how fluids interact with hairy surfaces, and here we aim to make progress in this direction. In 
particular, we ask how the elasticity, wetting angle and incline of an array of hairs determine their effectiveness in supporting a 
liquid away from the base substrate. 

A water drop deposited on a partially wetting solid does not spread indefinitely, but remains as a localised spherical cap, 
making a finite angle with the substrate. This angle is determined by the balance of surface tensions between water, air and solid 
surface, as given by Young's law Jfl 

a 7sa - 7sw .,s 

COS# Y = , (1) 
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where 7, 7sa and 7sw are the surface tensions of water-air, substrate-air and substrate-water contact. 8y, the Young angle, 
describes a hydrophilic, hydrophobic or neutrally-wetting surface depending on whether it is below, above or equal to 90° 
respectively. Using known materials, the Young angle that can be achieved for water on a smooth substrate is no more than 
120°. However some rough surfaces can produce much higher contact angles. This occurs when the drop lies on top of the 
topographical structure, in the suspended or Cassie-Baxter state 01. In the suspended state, the extra air-liquid interface at the 
base of the drop can result in a very high contact angle, even when the material is only slightly hydrophobic. Moreover contact 
angle hysteresis is small and the drop rolls over the surface very easily [1], so the suspended state is superhydrophobic. 

The drop can also penetrate the interstices between the relief, to lie in contact with all points of the substrate, in the collapsed 
or Wenzel state HI]. The collapsed state boosts the contact angle only modestly, and increases resistance to a drop's motion 0], 
so is usually not considered superhydrophobic. Ensuring the feasibility and robustness of the suspended state is therefore an 
important consideration in the design of superhydrophobic surfaces. 



Facilitating a suspended state for hydrophilic Young angles would also have important applications. A major weakness of 
many superhydrophobic surfaces is that they fail to repel liquids other than water, for example hydrocarbons. This is because 
there are virtually no materials for which 8y with respect to these liquids is greater than 90° QJ. Oily substances will sit in 
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the collapsed state or, worse, be imbibed into the structure. This will despoil the surface, because a water drop will now be in 
contact with the oil, rather than an air layer. If the surface can be made superoleophobic ifioll . sustaining a suspended state for 
Young angles below 90°, then this problem will be solved. 

Posts on a substrate with a degree of overhang may permit suspension, even when the Youn g angle is hydrophilic, by 
pinning the contact line IXltl - Such surfaces have been constructed using mushroom shaped posts ifloL 1 12fl . As an alternative 
to overhanging geometries, we investigate the possibility that flexible posts or hairs on a surface may deform to support water 
away from the surface, even when they are hydrophilic. 

The leaves of the Lady's Mantle (Alchemilla Mollis/Vulgaris) are superhydrophobic lfl3l [Till , with water forming beads 
on them. Closer inspection shows them to be covered with fine hairs, just visible to the naked eye. Surprisingly, experi- 
ments [13] found these hairs to be hydrophilic. Otten and Herminghaus 11311 stressed the importance of the flexibility of the 
hairs in sustaining the suspended state. They proposed that the hairs bunch together to reduce distortion of the interface, 
and the consequent energy cost of bending the hairs prevents them from being entirely wetted. These arguments were dis- 
puted by Berndino et al 1151 . who cast doubt on both the viability of the mechanism, and the elasticity of the Lady's Mantle hairs. 

Some arthropods, such as the water spider and the water boatman, can dive underwater for substantial periods. Observations 
show a shiny appearance to the underside of these animals flql. which is due to a layer of air trapped against the body, called 
a plastron. The plastron provides the means of respiration lfl7ll . not by acting as an artificial lung (a store of air taken from 
above water) but as an artificial gill (an interface through which gas exchange with the water may take place) fTill . Thorpe 
and Crisp II19I423I1 demonstrated that the air layer is trapped by a bed of hairs, such that the surrounding water is maintained 
in a suspended, superhydrophobic state. In contrast to the Lady's Mantle, it has been argued that flexibility undermines the 
water repellency of the plastron of a diving arthropod, with any bunching of hairs leaving gaps through which the water will 
invade i23il24ll . Conversely, it has been argued that flexibility of the hairs on an insect's legs can provide directional control as 
it walks on the pond surface l24ll . These examples suggest that hair flexibility modifies the wetting properties of a surface in a 
rich and nontrivial way. 

A mathematical model of the static configurations of a drop on a hairy surface can be constructed by treating the drop 
as a capillary surface and the hairs as elastica. In recent years, there have been a number of such elastocapillary analy- 
ses lfl5l l25l - [32ll . Of particular relevance to the work we present in this paper are those investigating the interaction of initally 
flat elastic sheets with an initially flat fluid interface. Neukirch et al 112811 . found that a perfectly wetting elastic beam, raised 
and lowered from a tank of liquid, exhibited hysteretic behaviour, arising from the metastable binding of the beam to the 
interface. Park and Kim [29] similarly showed that the threshold load for a clamped rod to pierce an interface was greater for 
a flexible rod than for a rigid rod. Kim and Mahadevan [30] studied the problem of capillary rise between a pair of elastic 
sheets arguing that, at modest flexibility, the height of capillary rise is increased by the narrowing of the gap, but for higher 
flexibilities the sheets can seal together, thus preventing further rise. These three studies l28l - [30ll take the capillary length to be 
comparable to the system size so that gravity is an important influence which sets a decay length for disturbances of the fluid 
interface. The study we present in this paper differs from these by considering a periodic system with a capillary length much 
longer than the periodicity. Gravity may therefore be ignored, leaving just a balance between the interfacial and bending energies. 

The model we treat is an interface supported by a line of flexible hairs. In general, the hairs will lie along the interface, 
thus alleviating the surface tension of the drop. We obtain an exact analytic solution in two dimensions which allows us to 
describe the possible drop configurations and obtain quantitative boundaries for the transitions between them. For vertical hairs, 
as for rigid posts, a drop placed in the suspended state can remain there only if the hairs are hydrophobic. However, for both 
hydrophobic and hydrophilic hairs, we identify a partially suspended state where the interface lies between the tops of the hairs 
and the surface. We show that, once the drop is in the partially suspended configuration it can remain there even for hydrophilic 
hairs. We discuss whether it is feasible for a drop placed gently on the surface to access the partially suspended state, and show 
that this is rather natural, if the hairs have a small angle of inclination to the vertical. 

The layout of this paper is as follows: In Sec. [II] we describe the model and characterise the possible stable and metastable 
interface configurations. In Sec. [Ill] we write down the free energy and obtain equations for the profile of the hairs. The stability 
of partially suspended singlet and doublet states, which have the periodicity or twice the periodicity of the lattice of hairs, are 
discussed, as a function of the hair flexibility and Young angle, in Sees. |IV| and [Vl respectively. We then, in Sec. [VI] present 
numerical results showing how the partially suspended states can be formed. Sec. I VIII considers inclined hairs, and Sec. I Villi 
discusses the effects of finite drop curvature. In Sec.[lX]we extend our results to three dimensions. A conclusion summarises the 
paper and discusses the many directions for further research. 



FIG. 1 : Illustration of the Pi state, showing some of the parameters and variables of the system. The hairs are shown as green (lighter grey) 
and the liquid-air interface as blue (darker grey). 
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FIG. 2: Schematic illustrations of the periodic states of the system. The hairs are shown as green (lighter grey) and the liquid-air interface 
as blue (darker grey). At the top is the fully suspended state, F, and at the bottom is the collapsed state, C. On the left are the singlets 
states: partially suspended, Pi, and quasi-suspended Qi. On the right are the doublet states: partially suspended, P2, quasi-suspended Q2, and 
unzipped U. See section|Il]for details. 

II. MODEL FOR A LINE OF HAIRS 

We consider an infinite row of hairs (or equivalently elastic sheets with no bending along y), labelled with integers n, attached 
to a flat substrate, with regular spacing D along the x direction, and with a fixed inclination il to the z direction (see Fig. [TJ). 
The hairs are inextensible, of length L and width negligible compared to D, and have a bending modulus K. They are in contact 
with the base of a liquid drop of surface tension 7. The surfaces of the hairs and the substrate are smooth and have a contact 
angle with respect to the liquid of the drop of 8y- The drop is assumed to be sufficiently small that gravity can be neglected, 
but large enough that Laplace curvature and edge effects are unimportant, (L, D) <C R <C Ac where R is the drop radius and 
Ac is the capillary length. In this regime the interface between neighbouring hairs is straight. Note that a hair may contact the 
interface in different places on its two sides. This is a consequence of working in two dimensions. We discuss the circumstances 
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under which our results are relevant to hairs in three dimensions in Sec.lIXI 

We introduce a classification scheme for the stable interface configurations, which are illustrated in Fig.|2] In the rigid limit 
it is well known that the possible states are the fully suspended or Cassie-Baxter state, in which the interface rests on the tips of 
the hairs, and the collapsed or Wenzel state in which the liquid wets the hairs and the base substrate. We denote these F and C 
respectively. The symmetry of the interface configurations implies that the hairs in either of these states will be straight for all 
values of the hair rigidity. We choose the free energy of the fully suspended state F to be zero. Therefore the free energy of the 
collapsed state, C is per length D along the surface. 

F c = -jD + ( 7sw - 7sa) (D + 2L) 
= -f[-D-cos8 Y (D + 2L)} 

where Eqn. (HJ has been used in the second step. 

Flexible hairs also permit partially suspended states, where the interface on one side of a hair is pinned at its tip, but that on 
the other side meets the hair part-way up, at an arc length sq from the base of the hair, say (see Fig. [TJ. The simplest partially 
suspended state, which we shall term the singlet and denote Pi, has periodicity D. Each hair bends in an identical curve, in 
the same direction, and a section of interface connects the tip of hair n to the point sq on hair n+ 1. A second possibility is a 
doublet, P2, with periodicity 2D, where each hair takes the same shape to within a reflection, but neighbouring hairs bend in 
opposite directions. In this state, the section of interface between hairs n and n + 1 joins the tips of neighbouring hairs for n 
even, and the points sq for n odd. Hence the interface always lies parallel to the substrate, but at two different heights. 

Note that F is equivalent to the partially suspended states for so = L. As so —> in the Pi configuration, the interface meets 
the substrate at an obtuse angle and the fluid readily wets the surface together with the dry sides of the hairs resulting in the 
collapsed state C. By contrast, as sq — > in the P2 configuration the substrate can wet between hairs n and n + 1, n odd, but 
remain suspended for n even, resulting in air pockets. We shall refer to such configurations as unzipped states, U (see Fig. |2). 
Other conceivable states are quasi-suspended states Qi,2 where the interface meets a hair part-way along on both sides. We 
shall show, in section [TV] that such states are not stable here. 



III. THE FREE ENERGY 



In this section we write down a free energy for the partially suspended states Pi 2. starting with the case O = 0. The shape of 
the hair is defined by the position vector r(s) where s is the arc length which runs from at the base to L at the tip. Because s is 
a measure of arclength, r(s) is the unit tangent to the hair at s and f(s) is the directed curvature. The free energy of the system, 
measured for a length D along the surface, and measured relative to the free energy of the fully suspended state, is 

T=\ [\Ky 2 + \<j{s)(y 2 -l)]ds + [\Ky 2 + \<j{s) (y 2 -l) - 7 cos 6> Y ] ds + 7 (|A| - D). (3) 

JO J s a 

The first two terms in each of the integrals correspond to modelling the hairs as inextensible Euler elastica. The first term denote 
a free energy quadratic in the curvature. (Such a functional form may be justified from the considering the hair to be made of a 
linear elastic material ll33ll . with K = EI, where E is the Young's modulus of the material and / is a geometrical moment of 
the hair's cross section.) The third term in the second integral is the free energy which results from the liquid wetting the hairs. 
The remaining term is the liquid-gas interfacial energy where 

f£>e x -(r(L)-r( So )) , for Pi , 
\[D- (x{L) - x(s ))} e x , for P 2 , 

is the length of interface between successive hairs. 

We split the minimisation of the free energy into two stages, first applying the calculus of variations for fixed so, and then 
minimising with respect to so. In the first step we consider the variation of the free energy T — > T + 8F with respect to 
infinitesimal perturbations of the curve r —> r + Sy. If the hair is in a configuration corresponding to an extremum of T then 
8F will be zero to first order in all Sy. However, we must constrain the set of solutions to those where the parameterisation 
corresponds to s being the arclength along the curve. This is important for two reasons: to ensure the hair has the correct length, 
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and to ensure that r corresponds to curvature. To do this, we also vary er(s), independently of r(s). a(s) may be likened to a 
stress in the hair, resisting extension or compression. Considering variation of T with respect to r(s) and a(s), 

f s ° f L A 5A 

8F = I {Kt St + ai.Si + (r 2 - l) 5a} ds+ { Kr.Sr + ar.Sr + (r 2 - l) 5a} ds + J-rrr- 
Jo Js„ \A\ 

= J°° | (KT - .5r + (f 2 - 1) fcj ds + [-K (r+ - f") .Si + (K (r + - r") - (a+ - a~) f +7) .6r] SQ 

+ J I (iCr - -St + (r 2 - l) 5a j ds + [Kr.5r + (-K'i' + ar - 7) .5r] L 

(5) 

where we have expressed the surface tension as a force vector 7 = 7jxj, and used the fact that r and r must be continuous if 
F is not to diverge. A continuity condition does not automatically apply to other quantities; hence the ± superscripts to denote 
the direction from which the limit s — > sq is taken. 

To find the extremal solution we set 5 J 7 = 0. Since 5r and 5a are considered independent and arbitrary functions, the 
coefficients of each of these in the integrand of (f5]l must be zero for all s for the condition 5 J 7 = to hold identically. From the 
coefficient of 5a, 

r 2 - 1 = , (6) 

or |r| = 1, thus guaranteeing that s measures arclength. a(s) may be viewed as a Lagrange multiplier which acts locally 
to constrain the curve derivative at each point (as opposed to a constant Lagrange multipler, which would impose a global 
constraint). 

Putting the coefficient of 5r in Eqn. <J3J to zero gives a differential equation for the shape of the hair 

K-t-«&=0, (7) 

ds 

with the boundary conditions 

r(0) = , (8) 

r(0) = e z , (9) 

Kr-(s ) = Kr+(s ) , (10) 

[K (r + - r") - (a+ - a~) v] sq = -7 , (11) 

Kv(L) = , (12) 

[-Kr +ai] L =i . (13) 

Eqns. d8l91 l are imposed boundary conditions, which arise from fixing the position and orientation of the hair at s = 
respectively. Eqns. ( 1 1 OH 1 3b are natural boundary conditions, arising from the boundary terms of Eqn. ©: 5r, 5i are free at 
s = so, L so their coefficients again must be zero. 

Integrating @, and using (fTTT i and (fT3l l to set the constants of integration, gives 

fo, < s < s , 
K r — ar = < (14) 
y -7, s < s < L . 

It is useful to note that the righthand side of the equation is the net force exerted on the hair above the point s. The unknown 
function a may be eliminated by dotting ( TBI with r and using r.r = ^d(r 2 )/ds = 0. Integrating a second time, and using 
Eqns. ( fT2b and ( fTUl ), gives 

iA^ 2 r^-/ T(i " (L)_i " (so)) 0<s<s °< fl « 

2 KV lSj -l 7 .(r(L)-r( S )) s < s < L . ^ 
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It is useful at this point to switch to an angular representation. Writing f = sin ip e x + cos ip e z and 7 = 7 (sin a e x + cos a e z ), 
and defining the elastocapillary length A = J—, Eqn. ( fT3T > becomes 



A 2 (cos [a — ^>(£)] — cos [a — ip(so)]) < s < s , 
A -2 (cos [a — V>(£)] — cos [a — ?/>(s)]) so < s < L 



(16) 



Eqn. dT6b implies that, in the lower segment of the hair, < s < so, VK 3 ) increases linearly with arc length 

if)(s) = X~ y2 (cos [a — ?/>(£)] — cos [a — ip(so)])s . 

For so < s < L the nonlinear differential equation for the shape of the hair can be solved in terms of elliptic functions |3 
Making the substitution 



(17) 



sinx(s) 



cos [\{a- ip(s))] 
cos 



in Eqn. ( TToT l gives 



2x 2 cos 2 x- 







cos 2 


'a-tp(L)' 
2 


sin 2 


a— ip(s) 
2 




x 2 



a-j>{L) 
2 



2A" 2 cos 2 



2 



a— i/j(s) 
2 



[1 - sin x) 



which may be integrated in terms of elliptic integrals 



A" 1 (L - s) 



= A ( 1 - cos 2 



dx 



a-ip(L) 
2 



sin 2 x 



1 — cos 2 



2 



. 2 

sin x 



^cos[i(a-i/>(L))] 



F cos[i( Q ^(L))] 



(18) 

(19) 
(20) 
(21) 

(22) 
(23) 



where K m and F m denote complete and incomplete elliptic integrals of the first kind with modulus m 114511 . The expression 
inverts as 



C11 cosri(a-'i/j(L))] 



[\-\L-s)] 



smx(s) 



dn 



cos[i(a-'i/j(L))j 

where cn m and dn m are the Jacobi elliptic cosine and delta. 



[X-HL-s)] 



(24) 



if>(s) must be continuous, so two implicit equations for the constants ip(L), ip{ s o) an£ l a follow from taking s = so in 
Eqns. ([171) and ([241 



ip(s ) 2 = 2 (y) (cos [a - V>(£)] - cos [a - ip(so)\) , 

[A-^L-s 

cos [§(a - i>(so))] = cos [5(0 - rp(L))] 



(a-^(L)) 



dn 



[A-i(i-«o)] 



(25) 
(26) 



The third equation determining these parameters follows from the geometry of the model. For the doublet state P2, symmetry 
implies a = ^, provided the hairs do not overlap. For the singlet state Pi the position of the contact points on successive hairs 
and the slope of the interface a are related by 



cot a = 



■ J^ o cos ip(s)ds 
D — smtp(s)ds 



(27) 
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which by trignometric rearrangment becomes 

D cos a+ / (sin [a — ^P(s)] cos 2 a — cos [a — V'(s)] cos a sin a) ds = — / (cos [a — ip(s)] cos a sin a + sin [a — ip(s)] sin 2 a) 

Js J s 

(28) 

Therefore, 

,L 

Dcosa — — / sin [a — ip(s)} ds , (29) 

J s 

which, by use of the first equality in Eqn. ( TToT l, may be written 

/^W Asin[q-V»(g)] 
Z?cosa = - / a-0 
•Ak«o) v 2 (cos [a - ip(L)} - cos [a - ^( S )D 

= -A-v/2 (cos [a - ?/>(£)] - cos [a - ^( So )]) OO) 

= A 2 ^(go) 
SO 

where we have used Eqn. ( |25T l in the final step. 

Substituting Eqn. ( [ToT l into Eqn. (0, integrating, and using Eqn. d24l >. the free energy may now be written. 

— =L cos [a — ip(L)] — so cos [a — tp(so)] — (L — s ) cos 9y + cos (a — ip(s)) ds + |A| — D 

7 Js 

=Lcos [a — ip(L)] — s cos [a — ip(so)] — (L — s )(cos9y + 2) (31) 

+ 4 ( E cos[i( a -^))] (5) - E c OS [I (Q ^ (i ))] W-o))) - U (1 - sma) , 

where E is the elliptic integral of the second kind. To find the minimum of the free energy for a given Young angle we now 
implement the second part of the extremisation process; taking the derivative of _F with respect to So- After some algebra we 
obtain 

dj~ 

- — = 7 (cos [a ~ tJj{s )] + cos 9y) , (32) 
ds 

which is zero when 

V^(so) - a + n = Y ■ (33) 

Hence, as expected, in equilibrium the angle between the interface and the hair is the Young angle. Calculating the second 
derivative of T leads to the expected condition for a stable equilibrium, that a small advance of the contact line must reduce the 
contact angle 

dQ Y 

-j— > • (34) 

ds 

Finally, we have assumed a partiallysuspended state which will only be stable if the interface remains pinned to the tip of the 
hair. Invoking the Gibbs' criterion 113511 . pinning will occur if 

a - tj){L) < 9 Y . (35) 



Thus the contact angle at the tip may take a finite range of values, as occurs, for example, on the corner of a grooved surface [3 
In conclusion, partially suspended interface configurations, where the interface runs between the tip of a hair and a point part 
way up the neighbouring hair, are found by substituting Eqn. d33l to eliminate ip(so) from Eqns. d25l ), d26| ) and ( f30l > 

(a + 9 Y - tt) 2 = 2 f ?f) 2 (cos [a - if)(L)] + cos Y ) , (36) 



A 



[X-^L-so) 



ri -i ri i cosf^a— i>(L))] 

sin [I0 Y ] = cos [l(a - m) ] dn * - s )] ' (3?) 

cos^(q— 0(L))J 

J tt — 9y — \~ 2 Ds cos a for Pi , 

(38) 



I 2 

These states are viable minima given that the conditions (f3~4-b and (l35T l hold. 



for P 2 
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FIG. 3: Variation of the reduced free energy .F/7 per length D, given by Eqn. J3U . with the position of the interface-hair contact so for the 
partially suspended singlet state Pi. Results for w = y^y/KL = 0.25 (red), l/y/2 (orange), 1.25 (green), 1.75 (blue) and 2.25 (magenta) 
are shown for different values of the equilibrium contact angle (a) 9y = 90°, (b) 8y = 110°, (c) 0y — 70° 



IV. SINGLETS 



We first assume that the drop is in the partially suspended singlet state Pi and discuss when this state can remain metastable 
with respect to the fully suspended and collapsed configurations F and C. A similar analysis for the partially suspended doublet 
is given in the next section. The results are summarised in Fig. |4(a)| and are presented in terms of the dimensionless parameter 
w := \~ 1 L which is a measure of the ability of the interface to bend the hairs; w = for rigid hairs. 

Consider first neutral wetting, 6\ = 90°. The variation of the free energy of the partially suspended state Pi is shown in 
Fig. |3(a)| for different values of w. Recalling that the suspended state F corresponds to sq = 1, and that its free energy is chosen 
as zero, the curves show that, for small w, F is stable. For w > wq, however, sq = 1 becomes unstable, and a new stable 
minimum, corresponding to Pi appears. Its position decreases continuously from so = 1 as w increases from wo signalling a 
continuous transition. In the Appendix we show that u>o = -J 'L/D. 

There is similar behaviour in the hydrophobic region of parameter space, as shown by the free energy curves in Fig. |3(b)| for 
6y = 110°. A drop in the partially suspended singlet state will rise to be fully suspended for small w, but will remain in Pi for 
larger w. The new feature for 8y > 90° is that the free energy has a maximum near sq = 1 indicating that the Pi — > F transition 
is not reversible. The boundary of metastability moves to higher iv with increasing 8y as shown in Fig. |4(a)| for D/L = 2. 

A typical free energy plot for hydrophilic hairs in the partially suspended state is shown in Fig. |3(c)| For hydrophilic hairs, 
Pi is susceptible to collapse by two distinct mechanisms, which place upper and lower limits of stability on the state. For large 
if there is no stable Pi state. The interface will descend to so = and then spread across the substrate to form a collapsed state. 
For smaller w the free energy curves have a minimum corresponding to a Pi state that is stable to collapse. The minimum moves 
continuously to sq = as w increases but the transition is not reversible as there is a strong free energy barrier associated with 
dewetting the base substrate. Alternatively, if w is decreased to 0, the minimum in the free energy OTb persists which is, at first 
sight, odd because Pj is not stable in the rigid limit. This is resolved by noting that the Gibbs' criterion (f35T > is violated for small 
w. The partially suspended singlets are no longer stable and there is a first order transition to the collapsed phase. Free energy 

3(c)| The explicit form of the depinning curve 



curves which correspond to these values of w are indicated by dotted lines in Fig 



may be found by substituting Eqn. d35l l (written as an equality) into Eqns. ( I36H38 



, and eliminating a and sq, 



L 1 



A 2Vcos6V 



-2^v/cos#y 



7\ 



K e Y — F Y (arcsinftan , (39) 

rne — — rne — — * Z / 



cos ~2~ cos 2 



where K rn and F m denote complete and incomplete elliptic integrals of the first kind with modulus m. The boundary where 
Pi becomes unstable to C changes direction at a prominent cusp marking the crossover between the two different collapse 
mechanisms. The cusp is significant because it marks the lowest (9y for which Pi is metastable for a given L/D. The position 
of the cusp can be found by substituting the depinning condition (l35T l and so = into Eqns. (1361138b . giving the equalities 

L 2L 

A = ^s^ - F cos^ (™ - JJ^- ^ 



We now describe how the metastability of the partially suspended singlet state varies with the hair spacing to length ratio D/L. 
wq, the position of the continuous Pi O F transition for neutral wetting, decreases as D increases thus making Pi states 
accessible to hairs of greater rigidity or, equivalently, to liquids of lower surface tension. In contrast, the range of hydrophilic 
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(a) 




(b) 

FIG. 4: Phase space diagrams for D/L — 2. (a) Black dashed lines show the limits of (meta)stability of the partially suspended singlet state 
Pi and red dashed lines the limit of stability of the fully suspended state F. Note that Pi can persist at hydrophilic contact angles, (b) Black 
dashed lines show the limits of (meta)stability of the partially suspended doublet state P2 and of the unzipped state U, and red dashed lines the 
limit of stability of F. Boundaries correspond to transitions which take place in the direction of the arrows, w = \Ay/ KL is a dimensionless 
parameter describing how easily the hairs can be bent by the fluid interface and 8y is the equilibrium contact angle of the hairs. Diagrams have 
been added which show the hair and interface profiles at given positions in parameter space: hairs are depicted in green (lighter grey) and the 
liquid interface in blue (darker grey). 
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Young angles for which Pi is metastable is reduced, with the cusp in the Pi — >• C curve moving to higher contact angles. 
Indeed, in the limit D — > oo, Eqns. JA71 > and d40l ) show that the cusp and critical point converge at (9y = 90°, w = 0), and there 
are no partially suspended states in the hydrophilic region. 

We consider another possible drop configuration, Qx, where both ends of a segment of interface lie part way down the posts. 
This state could, for example, be created when an interface in the Pi configuration depins from the tip of the hairs. In addition 
to the contact point at s = so, the second contact point in a Qx state will be situated at a higher point s = L. Above the two 
contact points the hair will experience no forces and hence will not be bent. Therefore equilibrium conditions for Qx can be 
obtained by making the substitutions L —> L in Eqns. d36ll38l ). However, the analogy is not complete because Eqn. ([35), which 
represented the Gibbs' criterion for P states, is now an equality. Thus Qx states obey Eqn. ( 1391 ). In order for these equilibria to 
be stable we require ^j- > 0, in analogy with Eqn. d34b . or equivalently that > 0. However, examination of the righthand 
side of Eqn. (l39l shows it to be a monotonically decreasing function of 9y for all values of X/D, and therefore there are no 
stable Qx states. 

Finally, having found stable equilibria for a systems in which periodicity is enforced, we should check that such states are not 
unstable to variations between hair profiles. To do this we consider a more general system, where each hair n has an independent 
profile, and check the stability of the uniform states, given by Eqs. ( I36H38I ), against perturbations of all possible wavelengths 
along the chain of hairs. We find stability for all Pi states, except for a small region along the F — >• Pi curve close to the critical 
point. For the system considered here, this region extends only ~ 1°. 



V. DOUBLETS 



We now analyse the stability of the partially suspended doublet states P2 against transitions to F, C and the unzipped state U 
(see Fig. [2j. These are easier to handle because the additional symmetry implies that the interface lies parallel to the surface, so 
a = w and the results are independent of D (provided that the hairs are sufficiently spaced not to touch). 

The equations of equilibrium follow from substituting a = ? into Eqns. (f36l > and d37| i, giving 



(0y - |) 2 = 2 (^) 2 (sin r/,(L) + cos 6 Y ) , (41) 



cn 

cos I 



cos^(-2-V(£))J 

A phase diagram for doublets, constructed by solving Eqns. (l4TT i and d42l . together with the stability condition ( [34-b and the 
depinning inequality d35l ) is shown in Fig. |4(b)| 



The hydrophobic portion of the phase diagram is similar to that for singlets: P2 undergoes an irreversible transition to F as w 
decreases. The region over which P2 is stable increases with decreasing contact angle. The partially suspended doublet exists 
only for hydrophobic hairs because sq — > as 9y — > 90°. At the 9y = 90° boundary P2 becomes unstable to the unzipped state 
U. This in turn collapses to C for smaller 9y, when the bridging sections of interface in U depin from the tips of the hairs. 



VI. PRODUCING THE PARTIALLY SUSPENDED STATES 

In Sees. [IV] and [V] we discussed the stabilty of the partially suspended states assuming that the system was initially in 
these configurations. We found that these states can remain stable against transitions to F, C or U over a substantial range 
of parameter space. In particular the partially suspended singlet persists at hydrophilic contact angles. We now address the 
question of how easily the partially suspended states can be created in the first place. 

Consider a large drop in the fully suspended state, F. If the Young angle is not hydrophobic, F is unstable and the system 
must undergo one of three available discontinuous transitions: to Pi, U or C. We find that, for vertical hairs, the chosen state 
depends very sensitively on any initial perturbation. We demonstrate this using simulations, where the fluid and hairs are 
modelled using the free energy lattice Boltzmann method and the lattice spring method respectively. Details of these techniques 
are given in IB7H39I1 . In Fig. |5(a)| the interface is initially placed in F for 9y = 90° and allowed to evolve with time. A random 
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(b) 

FIG. 5: Lattice Boltzmann simulations showing the final configuration of an interface which is gently placed in the suspended state F on (a) 
vertical hairs, (b) hairs inclined at an angle of 5° to the vertical. Even a small angle of inclination favours the formation of a uniform, partially 
suspended singlet state. Parameters are w = 0.98, 8y = 90°. 

sequence of regions of Pi and TU result as a consequence of machine noise. 

This suggests that to robustly produce the partially suspended singlet configuration it would be helpful to impose a geometry 
which favours the formation of this state. This can be done by tilting the base of the hairs with respect to the vertical. 

In Fig. |5(b)| we show results obtained by repeating the simulation in Fig. |5(a)| but for hairs initially tilted at 5°. The effect of 
even this slight inclination is to form a uniform partially suspended singlet configuration. We therefore now extend the analysis 
presented in Sees. |IV|and[Vlto show how the region of stability of the Pi configuration is changed by aligning the base of the 
hairs at an angle SI to the vertical. 



VII. INCLINED HAIRS 



We consider hairs which have a fixed inclination Q (in the positive x direction, say) from the vertical at their base. 
Note that the reflection symmetry that holds for = is now broken, and we must distinguish two sets of partially sus- 
pended singlets, Pi + and Pi_, in which the interface between hairs n and n+1 is pinned at the top of hairs n or n+1 respectively. 

For singlets a non-zero inclination angle can be accounted for by a simple modification of Eqns. (|36| | and d38l l 

(a + 9 Y - n + ^) 2 = 2\~ 2 sl (cos [tp(L) - a] + cos 6 Y ) , (43) 
a + 6 Y - 7r =F = -A~ 2 £>s cosa , (44) 

with (|37T i unchanged. The functional form of the free energy (l3TT l remains the same. 

Consider first rigid hairs, w = 0. Simple arguments based on the Gibbs' pinning criterion show that F is stable for 0y > 
90° + O, while collapse occurs for 9y < 90°. Between these values the free energy is locally minimised by Pi + with the 
unpinned end of the interface lying at 

cos \Q — 9y] 

s a = L- D \ YJ . (45) 

sin fry 

Pi_ states do not exist in the rigid limit. 

The stability of the partially suspended singlet configurations for flexible hairs is summarised by Fig. [6] for D/L = 2 and 
fi = 5°. Comparing Fig. |4(a)| for vertical hairs, an important difference is the stability of Pi + at iv = 0. For w < wq the 
transition between F and Pi + is continuous, and hence reversible, and occurs along the line 8y = 90° + tt. For w > wq the 
boundaries of metastability of the two phases are different: F becomes unstable to Pi+ along 8y = 90° + il whereas Pi+ 
remains stable against a transition to F for higher 8y as w increases. 

The hydrophilic region of the phase diagram is qualitatively similar for vertical and inclined hairs. However note that the 
Pi+ states are stable less far into the hydrophillic region than for the upright hairs. The Pi_ configurations are, unsurprisingly, 
nearly always of higher free energy than the Pi+. However they do extend, as metastable states, further into the hydrophilic 
region. 
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FIG. 6: Phase space diagrams for D/L = 2 and a hair inclination of 5°. Black and blue lines show the limits of (meta)stability of the partially 
suspended singlet states Pi+ and Pi_ respectively, and red lines the limit of stability of the fully suspended state F. Dotted (solid) lines show 
the boundaries for discontinuous (continuous) transitions, which take place in the directions of the arrows, w = y^-y/KL is a dimensionless 
parameter describing how easily the hairs can be bent by the fluid interface and 6y is the equilibrium contact angle of the hairs. Diagrams have 
been added which show the hair and interface profiles at given positions in parameter space: the hairs bending to the right are in the Pi+ state 
and those to the left in Pi_ . 




FIG. 7: Schematic illustration of the buckling tranisition when a Laplace pressure is applied across the interface, for an increasing w. (a) for 
w < wb the vertical hairs do not buckle, (b) for w > wb the hairs buckle but the contact line may remain pinned in an F state, (c) for a higher 
w, the contact line depins from the hair tip, leading e.g. to a P state. 

VIII. THE EFFECT OF LAPLACE PRESSURE 



Having determined the equilibrium morphologies for our system, it is interesting to consider how a pressure difference 
between the two phases affects the system. An excess pressure in the liquid phase could represent the hydrostatic pressure 
experienced by a diving insect, or result from a finite drop size: Laplace's law relates the curvature of the interface R^ 1 to the 
pressure difference across it, through p = ^Rr 1 . 

Departures of the interface from planar, together with the pressure directly applied to the wetted part of the hair, complicate the 
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FIG. 8: Surface evolver simulations showing the interface and hair profiles for a rectangular array of hairs with spacing D x — 2L, D y — 4a 
and L = 32a. w — 3 and 9y — 90°. The interface closely conforms to the bending profile of the hair, with little curvature along y. 

equations describing the hair profile. However two basic effects of Laplace pressure on the phase diagrams may be understood. 
Firstly, the transition line for F — > C, moves into the hydrophobic region [40[]. This is because the bulging interface meets the 
hair at an angle 90° + f3 (see Fig.[7](a)) where 

sin^ = , (46) 

so Gibbs' pinning is overcome for 6y < 90° + /3. For elastic hairs, the critical point is also shifted, to 90° + ft + £1, and the 
value of Wq is slightly changed. 



A second consequence is that F states will undergo buckling. This is because there is now a downward force pD acting on the 
tips of the hairs. The profile of the hairs in the fully suspended state obeys a differential equation analogous to ( [Tol l 



\K^ 2 {s) = pD (cos[tt - V>(X)] - cos[tt - tp(s)}) , 
which can be solved to give a condition analogous to (|26| | 

2 



(47) 



dn . 



[wy/2 sin/3] 



(48) 



When fi = 0, Eqn. (gUl implies that either tp(L) = or cn . r ^j(l) 1 [iy-\/2 sin pj — 0. The second condition is equivalent to 



;V2sin/3 = K sin j^j 



(49) 



which can be achieved when 



W > Wb 



V8sin/3 



(50) 



This is Euler's buckling formula. Above wb, the ip(L) = solution to Eqn. (l48b is unstable and the hairs in F will be bent. 
For the case SI 7^ 0, Eqn. (l48l i indicates that there will be a finite degree of bending for any finite w. Should the hairs bend 
sufficiently, such that ip(L) + f + /? > 9y, the interface will depin from the tip, leading to a F —> P transition. 



IX. A THREE DIMENSIONAL MODEL 



We now consider an infinite drop lying on a two-dimensional array of hairs. The hairs are anchored to the substrate at z = 
to form a rectangular lattice, with separation D x and D y in the x and y directions. It is no longer meaningful to model the hairs 
as having zero width, and they are ascribed a constant radius a. This three dimensional geometry includes additional physics 
that results from a connected, two-dimensional interface. However, to make analytic progress building on our previous results, 
we need to keep the artificial constraint that the hairs can bend only in the xz plane so that their spacing along y remains constant. 

Partially wetting cylinders are bound to an interface because it is energetically favourable for them to lie there rather than 
in either bulk phase. If a cylinder is placed on a flat interface then it will adjust its position relative to the interface to achieve 
the correct equilibrium contact angle 8y- Similarly elastic hairs can be bound to an interface, as the decrease in wetting and 
interfacial contributions to the free energy can be greater than the penalty due to bending. 

We restrict ourselves to geometries where the interface curvature along y can be neglected. In general, we expect this to 
be a good approximation for (A, L, D x ) large compared to D y as borne out by the simulated hair profiles shown in Fig. [8] 
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(b) 

FIG. 9: Phase space diagrams for hair spacing D x — 2L, D y — 4a and radius a = L/32. (a) Black dashed lines show the limits of 
(meta)stability of the partially suspended singlet state Pi and red dashed lines the limit of stability of the fully suspended state F. (b) Black 
dashed lines show the limits of (meta)stability of the partially suspended doublet state P2 and of the unzipped state U, and r ed dashed lines 
the limit of stability of F. Boundaries correspond to transitions which take place in the direction of the arrows, w = ^jDy/KL is a 
dimensionless parameter describing how easily the hairs can be bent by the fluid interface and 9y is the equilibrium contact angle of the hairs. 
Comparing with Fig.[4] we see that Pi and P2 can exist at much lower contact angles than for the two dimensional geometry, but only at higher 
w. 
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These results, obtained using Surface Evolver 14111 . show that for D x = 2L, D y = 4a and hair radius a = L/32, the interface 
conforms closely to the bending profile of the hair with little variation in the y-direction. 

Under this approximation, we can write down a free energy analogous to Eqn. ([3]) 



T = I [\Kv 2 + \o(s) (r 2 -l)]ds+ / [\Kr 2 + \a(s) (r 2 - l) - 7 cos 8y2a(it - 6y)] ds 



+ 7 



D y {\A\-D x )+f (D y - 2a sin 9y)ds 

J So 



(51) 



The parts of the free energy ( BTT l relating to bending are unchanged in their functional form. The wetting contribution, in the 
second integrand, takes account of the hair making an angle dy with the interface, such that an area 2a(n — dy) per unit length 
of the hair is wetted. The interfacial part of the free energy, the last term of dSTb . contains two contributions. The first is from the 
area of the interface between the tips of one row of hairs, and the contact point sq of the next row in the x direction, assuming 
negligible variation of the interface in the y direction. We recognise that this part is the same as the interfacial term in Eqn. (01, 
but with 7 replaced by -D y 7, reflecting the change of dimension. The second contribution arises from the interface in the gaps 
between neighbouring hairs in the y direction, again taking into account that the interface meets the hair at the angle dy. We 
note that this term may be assimilated into the wetting term, so that the free energy in Eqn. ( ISTI l may be rewritten as 



/ [\Kt 2 + \a{s)(r 2 -l)]ds + / [\Ky 2 + \a(s) (r 2 - l) - jD y cos C ] ds + jD y {\A\ - D x ) , (52) 

JO Jsn 



where 



cos (9c = — 1 + 2—— [sin6*Y + (ir — 6y) cos(?y] 

Dy 



(53) 



may be thought of as the Cassie angle of the wetted section of the hairs. Eqn. d52b is equivalent in form to the free energy of 
the two-dimensional system ||3}, but with the replacements 7 — > jDy (and thus w y/jD y /KL) and dy — > 8c- Hence 
the two-dimensional Eqns. d34|36ll37] i, with these replacements, may be used to describe the configuration of the hairs in the 
three-dimensional model. 



A different modification is, however, needed to the depinning condition ( f35T >. If the interface were to depin from the tips of 
the hairs and descend a distance — SL, then the hitherto dry part of the hairs would be wetted, and the interface between hairs 
would be destroyed, leading to a free energy change. 

5F = {jcos8 Y 2a8y + ~/(D v - 2asin<9 Y ) - 1 D y cos [a - tp(L)]} SL , (54) 

from which it may be seen that the condition for p inning of the interface on the tips of the hairs, 5J-/SL < 0, is given by ((35), 
but with 9y replaced by an 'antiCassie' angle 14211 

cos 8 C = 1 - 2-^- [sin 8y - 8 Y cos 8y] . (55) 

Dy 



Using this mapping, we show the regions of phase space where the singlet and doublet partially suspended states are stable or 
metastable in Fig [9] The results presented are for vertical hairs, Q, = 0. Comparing to the equivalent two-dimensional plot, 
Fig. |4] we see that the three dimensional system is much more successful than the two dimensional one in producing partially 
suspended states on hydrophilic materials. In particular, for doublets, the unzipping transition is shifted from 9y = 90° to 
8c = 90° corresponding, for the parameters used here, to 8y = 57°. Although the model predicts U states, we expect these to 
be unstable to C, as the hydrophilic base substrate will promote spreading around the hairs. A higher w is needed to stabilise 
the partially suspended states, and there is no critical point where the P x <R- F transition becomes continuous. 

As the interface bends to follow the profile of the hairs in the x direction, it must have opposite curvature along y to preserve 
the condition that the mean curvature of an equilibrium interface is zero. For small D y , this curvature does not induce significant 
distortion of the interface but, as the spacing between the hairs increases, the approximation of two-dimensional hair profiles 
will become inappropriate. Moreover, the assumption that the hairs are constrained to bend only in the x — z plane is artificial. 
Indeed, cylinders sitting on an interface are unstable to lateral perturbations l43ll . with distortions of the interface causing them 
to cluster. More work is needed to understand the interplay between the tendency to cluster, the curvature energy of the interface 
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and the elastic energy of the hairs. 



We now give estimates to relate the results to physical systems. The dimensions of the hairy surface of the Lady's Mantle are 
documented |4|], while the stiffness of the hairs is still not known. We use our theory to predict the elastic modulus that would 
allow the drop to be supported. We have noted that w ~ wo corresponds to the regime where the drop may be supported by 
hydrophilic hairs, giving 




Using measurements from the Lady's Mantle 0], of L ~ 1 mm, D x = D y ~ 0.5 /mi and a = 10 /xm, taking the surface 
tension 7 ~ 0.07 Nm _1 for water, and using the formula for K = \-kEcl a for a 113311 cylindrical rod, we derive the suitable 
Young's modulus to be 

1 D 2 L 

E ~ —. 100 MPa (57) 

a 4 

which is a feasible value. 



As a second example, we consider elastic posts, created articially by Mock et al |4|]. These had dimensions a ~ 100 /im, 
d ~ 500 fim and L ~ 500 /im, and were made of HEMA with E ~ 0.5 MPa liill . The posts were square in cross-section, 
hence K = ±Ea 4 Q. A gain taking 7 ~ 0.07 Nm , gives w ~ 0.5, which is comparable to wq ~ 1. Unfortunately Mock et 
al could not obtain conclusive results for the interface position because of swelling of the posts. 



X. DISCUSSION 



We have developed an elastocapillary model to quantify the wetting behaviour of a surface patterned with a regular array of 
elastic hairs. Our primary motivation was to identify superhydrophobic states asking, in particular, whether these could exist on 
hydrophilic hairs. 

We began by working in two dimensions, considering a sufficiently large drop that interface curvature can be neglected, 
supported by a line of hairs or, equivalently, an array of parallel elastic sheets. This geometry allows an analytical solution. 
We identified three different states where the interface remains suspended above the base substrate. In addition to the usual 
Cassie-Baxter configuration, stable for rigid hairs, two partially suspended states can occur, the singlet (Pi), where all hairs 
bend in the same direction, and the doublet (P2) where neighbouring hairs bend in opposite directions. The partially suspended 
singlet can remain (meta)stable for both hydrophobic and hydrophilic hairs, for the dimensionless parameter w = ^j/KL 
within a suitable range. Simulations showed that a drop placed gently on the surface is much more likely to fall into the state Pi 
if the hairs are rooted to the surface at a slight incline to the vertical. 

The analysis naturally extends to a three dimensional system of hairs, given the restrictions that the hairs bend only in 
one plane, and the curvature of the interface perpendicular to this plane is minimal. The spacing of the hairs in the third 
dimension results in the Young angle being replaced by a Cassie angle with the consequence that, in contrast to two dimensions, 
hydrophilic hairs can support both singlet and doublet partially suspended states. The mechanism that stabilises these states is 
the propensity of the hairs to bind to the interface together with their resistance to bending. This is in contrast to that proposed 
in 11311 . where the hairs pierce the surface. 

The calculations reported were checked using both lattice Boltzmann and surface evolver calculations and, in future work, 
we aim to use the simulation approaches to obtain a more complete picture of what happens in three dimensions. It is important 
to understand the effect of instabilities which may cause the hairs to cluster: even in the rather artificial geometry considered in 
Sec. [IX] simulations show a tendency for hairs to bunch together in the y direction, and for more general hair configurations, we 
might expect more complicated clustering behaviour. Even the question of whether clustering bolsters lfl3ll or hinders ll23l l24ll 
superhydrophobicity still needs to be resolved. 

Another question to consider is the edge effects arising from a finite size drop. For example, in the case of Lady's Mantle, the 
hair spacing can be of the order of the drop size. We speculate that hairs will most easily lie in the interface at the edges of the 
drop, and that this may help the system to form partially suspended states. It will also be interesting to investigate the dynamics 
of drops on hairy surfaces to ask, for example, whether the flexibility increases or decreases resistance to motion, and whether 
singlet states may provide a switchable anisotropy of the surface. 
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Appendix A: Critical point for the partially suspended singlet at neutral wetting 

We already know that, for O = 0, the critical point lies on 9 Y = 90°, and we now derive an analytical expression for the 
critical value of w. 8 Y = 90°, so = 1 is always a solution to the equations for equilibrium (!36H38b . By differentiating these 
equations with respect to the contact parameter so, we can determine whether the equilibrium is stable or unstable, using the 
criterion (|34| i. The derivatives of Eqns. (|37| i and ((38), are 



d6 Y 
ds 

d0 Y 
ds 

d cn m x 



dip(L) 



s =L 



s =L 

= and 



ds 

da 
ds 



s =L 



da 
ds 



Sn=L 



(Al) 
(A2) 



sa=L 
d cn m x 



In deriving JA2b we use the results that ^ 

° * ' nx. dn T)1 x am da m x 

the remaining equation of equilibrium (1361 1 cancel at sq = L, so the second derivative is taken to give 



at x = (see e.g. 14510 . The terms of the derivative of 



dtp(L) 



ds 



sq=L 



2 r2 



ds 2 



d 2 ^{L) 



so=L 



ds 2 



To eliminate ^-pf and ^-J, we take the second derivative of (l37l i giving 



/ d9 Y 




da 


\ ds 




ds Q 


s =L 



d 2 Y 



ds 2 



Solving ( IA3HA4b yields 



Thus, the condition for stability (f34T > becomes 



2A" 



d8 Y (L) 



S =L 



ds a 



d9 Y 
ds 





(Pa 




S =i 


dsl 





(A3) 



s =L 



ds 2 



s =L 



d 2 a 
ds 2 



(A4) 



s =L/ 



D~ 



\~ 2 L. 



(A5) 



s =L 



A < A := VLD 



(A6) 



or, in terms of the dimensionless parameter w, 



W > Wo 



(A7) 
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